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1 Introduction 



Experiments with measurements of bremsstrahlung (Br) spectra in the a-decay of the nuclei 210 Po, 214 Po, 226 Ra 

and 244 Cm have caused an increased interest. One of the key ideas of the fulfillment of such experiments consists in 
finding a method of extraction of a new information about a-decay dynamics from the Br spectra (and a detailed 
information about the dynamics of tunneling). One can note a certain difference between the Br spectra pQ and 
03 for 210 Po, obtained experimentally for the values 90° and 25° of the angle between the directions of the a-particle 
propagation and the photon emission (these experiments and the difference of their spectra are discussed in @E]). 
One can explain such difference between the Br angular spectra on the basis of the following idea: the Br intensity 
depends on the directions of emission of the photons and motion (with possible tunneling) of the a-particle relatively 
the daughter nucleus. In such a way, a three-dimensional picture of the a-decay with the accompanying Br in the 
spatial region of nuclear boundaries has been devised. 

However, if the Br intensity varies enough visibly with changing the angle value, then one can suppose, that 
the photon emission is able to influence essentially on the a-decay dynamics and, therefore, to change all of its 
characteristic. From this point of view, the discussions open a way for obtaining a new information about the 
a-decay — through the angular analysis of the Br during the a-decay. But for such researches a model describing the 
Br in the a-decay, which takes into account the value of the angle between the directions of the a-particle propagation 
(or tunneling) and the photon emission, is needed. 

In theoretical aspect, some progress has been made here. One can note models of calculations of the Br spectra in 
the a-decay, developed on the basis of quantum electrodynamics with use of perturbation theory: the first paper [H] 
where a general quantum-mechanical formalism of the calculation of the Br spectra in the a-decay is proposed and 
the Br spectrum for 210 Po inside the photons energy region up to 200 keV was estimated (even until the fulfillment of 
the first experiments) ; essentially improved models in the dipole approximation [711$] and in the multipolar expansion 
[§] of photons current (wave function) with application of the Fermi golden rule; an approach |10j of the calculation 
of the Br spectra with realistic barriers of the a-decay) , models [TTlll2lll3 |l%] developed in semiclassical approximation 
(see also the Br spectra calculations in jH]), instant accelerated models [T]|^] constructed on the basis of classical 
electrodynamics (see also ^2]), methods [T31|14II15II11II12| . directed on a nonstationary description of the a-decay with 
the accompanying Br and the calculations of such non-stationary characteristics as tunneling time. One can recall 
also papers |lfiU17lll8lll9| with study of dynamics of subbarrier tunneling in the a-decay; an effect, opened in |2T)] 
and named Miinchhausen effect, which increases the barrier penetrability due to charged-particle emission during its 
tunneling and which can be extremely interesting for further study of the photon bremsstrahlung during subbarier 
tunneling in the a-decay). However, one needs to say that at this stage the calculations of the Br spectra by all these 
approaches are reduced to obtaining their integral (or averaged by angles) values and, therefore, they do not allow to 
fulfill an angular analysis of the experimental Br spectra (here, one can quote an approach in |14j based on classical 
electrodynamics, which shows a way for obtaining the angular spectra (see (15) and (17), p. 999), however here we 
shall use the direct quantum-mechanical approach of the Br spectra calculation, which describes the quantum effect 
of the subbarrier Br more precisely) . 

In [2] we had developed a multipolar method which takes into account the angle between the directions of the 
a-particle propagation and the photon emission. But the angular integrals used in this method are difficult enough to 
be obtained and some approximations are used, there was a convergence problem with the calculations of multipoles 
of larger orders, only the angular dependences of matrix elements of the selected multipoles El and Ml were found, 
while it is interesting to know also the angular dependence of the total Br spectra. Moreover, computer calculations 
of the Br spectra and their angular analysis will be essentially more complicated, if one passes from the type of the 
potential used in this paper (and also in 012]) to realistic potentials. In this sense, the approach proposed in [21] is 
not convenient enough. 

In this paper we present a new method for the angular calculations of the Br spectra in the a-decay of nuclei (started 
in |22| . with the resolution of a convergence problem in the Br spectra calculations existing in |22jV In our approach 
we introduce a simplified transformation, which reduce the complicated angular formalism of the calculation of the 
Br spectra (presented in [22) to a maximally simple form (this makes the method clearer and more comprehensible), 
with keeping the calculating accuracy as good as possible, where only one angle from all angular parameters is used 
— the angle used in experiments |3][5] for 210 Po. We show, that this proposed transformation works in the low-energy 
region of photons and, therefore, it can be applied to the analysis of all up to now existing experimental data of Br 
spectra in the a-decay of spherical nuclei (we relate 210 Po to them). In this paper we present the results of calculations 
of the angular spectra of the Br in the a-decay for 210 Po (the Br angular spectra in the a-decay have been obtained 
for the first time). We show, how the Br spectra are changed after deformation of the form of the a-decay barrier as a 
result of the correction of a component of the a-nucleus potential of nuclear forces (we did not find such calculations in 
other papers). An analysis of the problem of the calculations convergence of the Br spectra in the a-decay is included 
in the paper, whose resolution plays a key role for obtaining of the reliable values of the spectra. 
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2 A formalism of the calculations of the bremsstrahlung spectra in the stationary approach 

We shall consider the decay of the nucleus as the decay of the compound quantum system: a-particle and daughter 
nucleus. The a-particle is the electrically charged particle and during its motion inside the electromagnetic field of the 
daughter nucleus it emits photons. The spontaneous emission of the photon changes a state of the compound system, 
which is described by its wave function. For a quantitative estimation of the photon Br we use a transition of the 
system from its state before the photon emission (we name such state as the initial i-state) into its state after the 
photon emission (we name such state as the final f -state). One can define a matrix element of such transition of the 
system and on its basis find the Br probability during the a-decay (for convenience, we denote it as W(w)). According 
to [21], we obtain: 

W(w) = N k f w\p(w)\\ N = ^ cS l , k i)f = y/2mE iif , w = £k - E f , (1) 

where p(w) has a form 

p{w) = e(a) * f dr f rV/( r K ikr -|-^(i-)<OT. (2) 

a=l,2 { r 

Here Z e g and m are effective charge and reduced mass of the system, Eij, kij and ipij(r) are total energy, wave 
vector and wave function of the system in the initial i-state or in the final /-state (in dependence on the index i or / 
in use), ipi{ Y ) an d V'/M are the wave function of the system in the initial i- and the final /-states, e( Q ) is polarization 
vector of the photon emitted, k is photon wave vector, w = k = |k| is photon frequency (energy). The vector e'"' 
is perpendicular to k in Coulomb calibration. We use such system of units: h~ = 1 and c = 1. Notations are used in 
accordance with j2J- Similar expressions for the Br probability are used in |7||^H^) with further application of Fermi 
golden rule. 

In accordance with main statements of quantum mechanics, the wave functions of the system in the states before 
and after the photon emission are defined inside all space region, including the region of the subbarrier tunneling. A 
definition of a matrix element of the transition of this system requires an account of all space region of the definition 
of the wave functions of this system in two states. Therefore, we should include the tunneling region into the definition 
of the matrix element of the Br, irrespective of we know, whether the photons emission is possible during tunneling or 
not. 

Let's consider a subintegral expression in Here, the wave function tpi( T ) f° r the initial i-state and the wave 
function tpf(r) for the /-final state take into account the directions of propagation (or tunneling) of the a-particle before 
the photon emission and after it, correspondingly; the photons wave function (its main part consists in the exponent 
exp (— ikr)) points to the direction of propagation of the photon emitted. We see, that the quantum mechanical 
approach for calculation of the Br spectra initially has a detailed angular information about the process of a-decay 
with the accompanying Br. 

However, we see that a further development of the approach for the calculations of the Br spectra in the a-decay 
on the basis of the formulas QJ and @ by other authors (which consist in the calculations of p(w)) gives rise to 
angular averaging of the spectra. And a necesity has been arising in construction of an approach, which allows simply 
enough to calculate the Br spectra (with possible resolving the convergence problem in the computer calculations) 
with taking into account of the angle between the directions of the a-particle propagation (with possible tunneling) 
and the photon emission and without (essential) decreasing of the accuracy. 



3 A simplified angular method of the calculations of the matrix element 

In [22] the approach for the calculation of the Br spectra, allowing to find a dependence of the total Br spectra on 
the angle between the directions of the a-particle propagation (or tunneling) and the photon emission, was proposed. 
However, further research has shown, that it is extremely difficult to achieve a convergence in the computer calculations 
of the Br spectra by such approach and, therefore, such a method requires an essential development. Here, we propose 
a consecutive statement of such approach with a resolution of the convergence problem. 

Let's rewrite vectors e a of polarization through the vectors £_i and £ + i of circular polarization with opposite 
directions of rotation (see [2H], P- 42): 

e-i = 4( el - * e2 )> £n = + * e2 )- ( 3 ) 
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Substituting these values into J2J), we obtain: 



ft=— 1,1 

where 



+00 

p(w) = M2 / dr J r^r f (r)e- ikr ^Mr)dn, (4) 



h-x = 4=(1 - «), = — t=(1 + i), h-i +hi = -iy/2. (5) 
v 2 v 2 

Using the following properties (see [23i P- 44-46, |21|): 

^,W^^T 01 ,K), 
or dr 

T 01>0 (4)= 5j (1101-/^0)^,^(11;)^, (6) 

(U=-l,l ^ ' 

(110|1, -1,0) = (110| - 1,1,0) = Ji, 

where ( 1 10 1 — /i/iO) are Clebsch-Gordan coefficients and Tjj/ ^(n) are vector spherical harmonics (see [23], P- 45 and 
we use quantum numbers i = m = in the initial i-state), Yi jAl (n^) are normalized spherical functions (see 0], 
p. 118-121 (28.7), p. 752-755), we obtain: 

m'=-i.i 

Taking into account Q, and the ortogonality condition of the vectors and £ T i of polarization, we find: 

— +00 

K«) = -y~ E ^J^/W^/^^i^-^w)^ ( g ) 

where V/(r) = V/(^)^',m'( n r)- 

Let's consider the vectors k and r. The vector k is an impulse of the photon, pointed out the direction of its 
propagation. The vector r is a radius- vector, pointed out a position of the a-particle relatively a center of mass of 
the daughter nucleus and (because of mass of the daughter nucleus is larger sufficiently than mass of the a-particle) 
pointed out the direction of its motion (or tunneling). Then an angle between the vectors k and r (let's denote it as 
(3) is the angle between the direction n r = r/r of motion (or tunneling) of the a-particle and the direction n p h = k/fc 
of a propagation of the photon emited, i. e. it is the angle used in the experiments [3 | 3 f l]. One can write 

cxp(— ikr) = exp(— ikr cos (3), k= |k|, r = |r|. (9) 

Now we make such assumptions: 

— the photon emission process does not change the direction of motion (or tunneling) of the a-particle: 

nj, = n f r , (10) 

— the angle (3 is not depended on the direction of the outgoing a-particle motion from the nucleus region. 

Then, taking into account these assumptions and the ortogonality property of the functions Y; m (n r ), we obtain the 
following expression for p(w, (3): 

+00 

P {w,f3) = -sJ\ £ hj r 2 r f (r)^e- ik ^dr (11) 

and selection rules for quantum numbers I and m of the final /-state: 

the initial state: U = 0, rrij = 0; , 10 s 
the final state: If = 1, m/ = —fi = ±1. ^ ' 
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4 Spherical wave expansion 

For further computer calculations of the integral let's use an expansion of the plane wave in the spherical waves 
(for example, see j2U, P- 144, (34.1)): 

^.^ m+1W „0,(l)'(i')'^, (13, 



where z = rcos/3. Introducing spherical Bessel functions (see [21], p. 139, (33.9), (33.10) and (33.11)): 



/ 1 d \ l sin kr 



we obtain: 



and from wc find: 



^— ikr cos j3 I ^ikr cos j3 



= £7(-r/(2Z + l)fl(cos/3)j,(fcr) (15) 



r-+oo 

p(u;,/?) = -J-^i i (-l) i (2; + l)P i (cos/3) J] h„J mf (l,w), (16) 



where 



3 

2=0 M=-l,l 



•/».,(/•«•)= / r 2 ^)^^i*(^)dr. (17) 



J mf (l,w) is a radial integral, not depended on the angle /3. Now we obtain an explicit analytical dependence of the 
matrix element p(w , /?) on f/ie tmgfe /3 between the directions of the a-particle propagation and the photon emission (for 
the first time, in |21| the angular integrals were obtained for selected components of this matrix element — multipoles 
El and Ml only, with use of more difficult calculations). 



5 The bremsstrahlung in the Coulomb field 

Practically, in the numerical calculation of the Br spectra it is convenient to divide the whole region of the integration 
into two parts: the region 1 of a joint action of the Coulomb and nuclear forces not far from the nucleus and the region 
2, in which one can neglect by the action of the nuclear forces in comparison with the action of the Coulomb forces. 
Our analysis has shown, that an attainment of the convergence of the Br spectra calculations (which determines their 
accuracy, reliability of the found Br spectra) is reached first of all by correctness of the calculations in the region 
2. Namely, in this region it needs to solve a problem with definition of the external boundary of integration (its 
increasing leads to increasing of the accuracy of the obtained spectra, but to increasing of difficulty of the calculations 
and analysis), to choose the most effective method of the numerical integration (of an improper integral with an 
oscillating and weakly damping sub-integral function), to solve a problem with attainment of needed accuracy and 
convergence of the calculations. It defines time of the calculations, minimization of which appears extremely important 
for fulfillment of the real analysis of the obtained Br spectra in dependence on needed parameters. Therefore, maximal 
simplification of the formulas for the Br spectra in the region 2 is useful. 

Let's assume, that the potential, used in the radial integral 1|17JI . in spatial region of r is Coulomb since the value 
R c . We accept R c as the internal boundary of the region 2. One can write the radial integral J(l,w) in (|17|l so: 

J mf {I, W) = J in .rn s {I, w) + J C {1, w), (18) 

where 

Jin,m f (l,w) = J r 2 ip* f (r, to/) ji (kr)dr, 

°+oo (19) 

J c (l,w) = / r 2 ip* f (r) 1 ji(kr)dr. 

Ra 
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The radial integral J c (l, w) does not depend on the quantum number to of the systems in the final /-state. Then, one 
can write p(w,/3) so (with taking into account (JjJJ for the Coulomb component): 

p(w,0)=p in (w,P)+p c (w,0), (20) 

where 

■ +00 



Pin (w,j3) =-J^^i\-l) l (2l + l)P((cOS/3) J2 h ^n,m f (l,w) : 
' 2=0 ju=-l,l 

[2 + °° 

Pc(w,P) = W-^ i ' +1 (-l)'(2; + l)P ; (cos/3)J c (/,«;). 
' 1=0 



(21) 



We see, that there is no any interference between the components pi n (w,(3) and p c (w,(3) in the calculations of the 
total value of p(w,/3), but it exists in calculations of the total Br spectra. 



Then at I = we find: 



6 The first approximation at I = 

Legandre's polynomial of the order I equals (for example, see [21], p. 752 (c.l)): 

Pi(0) = ^J^yi 02 ~ !)'» MO) = 1, Pl(0) =0,0 = cos/?. (22) 

y m=-m (23) 

If for nuclei 210 Po, 214 Po, 226 Ra to use the potential with parameters as in [2] (and as in |7||21 also), then we 
find, that Br from the internal spatial region till to R c is extremely small (p in (w,f3) << p c (w,(3)). According to our 
estimations, for such potential it is smaller in 10~ 22 -10~ 24 times then Br from the external region. Therefore, one can 
neglect by Br from the internal region, and the total Br can be determined by the Coulomb field inside the barrier 
region and the external region. From (TJ3) we write down the Br probability in the first approximation at I = 0: 

W l=0 (w) = N k f w\pi l =°\w,l3)\ 2 = ^iV fc /U ;|j c (0,u;)| 2 . (24) 

One can conclude (it has obtained for the first time): 

— The Br probability in the first approximation at I = 0, formed by the Coulomb field both with taking into account 
of the nuclear forces of any shape, and without such forces, does not depend on a value of the angle (3 between the 
directions of the a-particle propagation (or its tunneling) and the photon emission. 

— The Coulomb field is degenerated by the quantum number m. This property distinguishes the Coulomb field from 
the nuclear forces at their account in the model. This difference is shown in the matrix elements (|23|) . The nuclear 
forces participate in formation of the decay barrier and, therefore, one can consider approximately them as forces 
working in the spatial region of the barrier, where there is a tunneling. One can assume, that one can divide the 
emissions from the barrier region and from the external region on the basis of the quantum number to/ . It can be 
interesting to find a possible way of extraction of the Br spectrum from the barrier region (or from the external 
region) from experimental Br spectrum on the basis of this property. 

7 The second approximation at / = 1 

Taking into account l|22l) . we find the Br probability in the second approximation at I = 1: 

pSi J) (tD,/5) = l\/3c0S/3 S~] /lpJm. m/ (l, W), 

m=-i,i (25) 
p { c =1 \w,P) = v^cos /3 J c (l,w). 
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Neglecting by Br from the internal region, we obtain the following expressions for the component of the Br 
probability of the second approximation at I = 1 : 



W (l=1) {w,(3) = N k f w p ( c l=1) (w,(3) = 6N k f w J c (l,w) cos 2 /3 
and for the total Br probability in the second approximation at I = 1: 

,J C (1,U>) 



Wi = i{w,(3) = Wi =0 {w) 1- N{w) cos P , N(w) =3i 



Jc(0,w) 



(26) 



(27) 



One can conclude (it has found for the first time): 

— The dependence of the Br probability in the a-decay on the value of the angle (3 between the directions of the 
a-particle propagation (or its tunneling) and the photon emission has harmonic type J2ZJ- 

— The account of the nuclear forces does not change the dependence of the Br probability in the second approximation 
at I = 1 on such angle value. 

— Exp. Ij27(l allows to find analytically maximums and minimums in the Br spectra in dependence on the angle (3. 



8 Convergence of calculations in asymptotic region 

There is an essential difficulty in the calculations of the Br spectra for the given nucleus, concerned with obtaining 
of the radial integrals (|19fl (or (|17fl ). This difficulty is caused by that such integral is improper, and its sub-integral 
function is oscillated and damped slowly with increasing of r. The function damps weaker with increasing of r, the 
larger region of integration should be taken into account in the numerical integration. For 210 Po the damping degree 
of the sub-integral function is such as for reliable values of the first 2-3 digits for the Br spectrum it needs to take into 
account (with the higher accuracy of calculations) 1 million of oscillations of this function. 
As an evident demonstration of this problem, let's consider an one-dimensional integral: 

+ OG 

sin a; , 

dx. (28 

x 

An exact analytical value of this integral at a = is known from theory of functions of complex variables, equal to it/2. 
The numerical calculation of the integral (with use of simple method of trapeziums, method of Gauss or other methods 
of the numerical integration) allows to obtain quickly the same result also, but with a given degree of accuracy (which 
determines a region of the numerical integration) . It proves a convergence of the computer calculations of such integral 
with a concrete choice of the parameter a. But at weak increasing of the parameter a the region of the numerical 
integration for obtaining of the same calculating accuracy for the integral $28\) increases essentially, and, therefore, a 
difficulty to calculate this integral numerically increases essentially. However, the application of methods of theory of 
functions of complex variables makes the calculation of such integral as simple again. So, on the example of the simple 
integral (|28|l one can meet with the numerical problem of the convergence of the calculations of the improper integrals 
with the damping slowly, oscillating sub-integral functions. 

We fulfill an analysis of the convergence of the calculation of the integral (|19fl on the basis of the analysis of 
damping of its sub-integral function in the asymptotic region, which is defined by wave functions in the initial i-, final 
/-states and the spherical Bessel function of order /. 

For enough large values of r one can use an asymptotic representation of the spherical Bessel function of order I: 

J i r ) (kr) = -t^(kr-^ (29) 

or 

An\kr) =(-l)« J ^(fcr) = (-ir Si " /,r 



, (30) 

where n is a natural number. 

The wave functions ipi( r ) an d V'/( r ) °f the initial i- and the final /-states are linear combinations of the Coulomb 
functions Fi(r),p) and Gi(r/,p) (divided on Pij, with quantum number I — or / = 1 for the initial i- or the final 
/-state, correspondingly). One can write the Coulomb functions for / in the asymptotic region so: 

F l {r],p)=sm0 l , Gi(>q,p) =cos0j. (31) 
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where 

Oi = P - V log 2/0+ + 07(77), p t j = hjr. 



2 , " mv (32) 

07(77) = argr(i77 + Z + 1), r}ij = - — . 

where .T(x) is Gamma function with argument x, v is Zomcrficld parameter. 
Now one can conclude: 

— The spherical Bessel function j^ as ^ (kr) in the asymptotic region damps (and oscillates) with increasing of r equally 
for any order I. 

— The Coulomb functions Fo(j]i,pi) and Go(rji,pi) of order for the initial z-state and the Coulomb functions 
Fi(rjf,pf) and Gi(rjf,pf) of order I = 1 for the final /-state damp in the asymptotic region with increasing of r 
equally, oscillate equally and are shifted at a phase between each other. 

— The total sub-integral function of the integral (|19fl in the asymptotic region damps with increasing of r equally for 
any order I. 

Taking into account (| 1 and we obtain: 

J { c as) (2n,w) = (~l)"ji as) (0, W ), J c (as) (2n+ l,w) = (-l)"J c (as) (l, w). (33) 

I. e. one can reduce any integral inside the asymptotic region to one of two integrals j( as ) (0, w) or j( as ) (1 , w) . 
Let's find the matrix element pi as {w,0) in the asymptotic region: 



pi^ (w, 0) = i 



\j\^ as) (0, u;) ^£(4n + l)P 2n (cos/3)^ + 
+ ^A a8) (l,w) ^(4n + 3)P 2 „ +1 (cos/3)^ . 



(34) 



Thus, we reduce the formula (jl 1|> for the Br spectra to the linear combination of two radial integrals, which are 
convergent (one can calculate them with a desirable accuracy limited by the calculations accuracy of a concrete 
computer) and do not depend on the angle, and factors - sums on n, into which the problem of convergence is carried 
out (one can meet with it in jlljl). 

It should seem, that one can cut off the region of the numerical integration in one boundary R for calculation of 
the integral J c (l,w) from (|19|l for any I. However, the calculation convergence of the integral is determined not only 
by the damping of the sub-integral function at large r, but by its behavior on the whole integration region also. An 
analysis has shown, that the sub-integral function inside the barrier region and inside the external region closer to the 
barrier behaves so, that the calculation of the total integral becomes more and more sensible to it with increasing of 
I and the calculation convergence becomes worse. Therefore, for obtaining the reliable values of the integrals J c (l,w) 
(for the same accuracy) it needs to increase the external boundary R of the integration region for larger I (such a 
conclusion has obtained by us is calculating of the angular Br spectra for 210 Po also). 



9 Angular calculations for the Br spectra in the a-decay of 210 Po 

As a demonstration of the described above method, let's calculate the angular Br spectra in the a-decay of 210 Po. 
For a comparison of results obtained in such approach, with results obtained by models |7||^||B). we shall choose the 
potential parameters as in |21 j (they coincide with the parameters of the potential with the external Coulomb field 
in and in 0). 

In spite of the fact that there are methods allowed to calculate absolute values of the Br spectra, in this paper at 
first we shall find the relative values of the Br spectrum for the given nucleus and then we shall normalize obtained 
spectra at one selected point of the Br experimental spectrum for the given angle value. This approach as against 
the previous one allows with a larger accuracy to analyze a behavior of the Br spectra in dependence on the angle 
(besides, it is more easy in application). 

In the beginning we calculate the total Br probability in the second approximation at I — 1 for the angle 90° by 
(|24|) with taking into account l|24|l (because the component of the Br probability (|26|l in the second approximation at 
I = 1 equals to zero at such angle) . Then we normalize the obtained spectrum by the third point of the experimental 
data ^ (we have such values w — 0.179 keV and W = 10.1 • 10~ 10 1 / kcV / decay), which were obtained for the 
angle f3 = 90° also. Knowing the normalized factor and using formulas (|27|) . we find the Br probability in the second 
approximation for the other values of the angle /3. 
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Table 1. Angular values of Br probability in a-decay of 210 Po in approximation I = 1 



The angular values of the Br probability in the second approximation at Z = 1 are shown in the Table 1. Here, 
one can see a variation of the Br probability in dependence on the angle /3, however this change is extremely small. 
The Br probability in the first approximation at / = coincides with the Br probability in the second approximation 
at I = 1 for the angle 90°. One can see that a contribution of the Br probability in the first approximation into the 
total spectrum is the largest for any angle value, i. e. it is extremely larger then the contribution of the component of 
the Br probability in the second approximation at I — 1 into the total spectrum. This conclusion has a physical sense 
(obtained for the first time): Br in the a-decay for 210 Po depends extremely weakly on the value of the angle between 
directions of the a-particle propagation (or tunneling) and the photon emission (in the given approach). The account 
of not zero component of the Br probability in the second approximation (for the angle values which are distinct from 
90°) increases the total Br probability. Absolute and relative variations of the Br probability relatively its maximal 
and minimal values 

AWi{w) = Wi^ x {w,(3 = 0°) - W l=1 {w,(3 = 90°); 

{w = |Wi =1 (n;,/3 = 0°)-Wi =1 ( W ,/3 = 90°)| (35) 

are included into the table also. 

The results of the calculations of the Br probability in the a-decay of 210 Po in the second approximation at I = 1 
by our approach are shown in Fig. ^ From the figure one can see, that for the angle 90° our method gives the Br 
spectrum, which is very close to the Br spectra, obtained by the models and However, as against the models 
[7|0|| our method shows the angular variation of the Br spectra. 



10 An analyses of the maximums and minimums in the angular Br spectra 

Let's find the values of the angle (3 between the directions of the a-particle propagation (or tunneling) and the photon 
emission, at which the Br probability has the maximal and minimal values. Using a derivative 



dW l=1 (w,0) 



Wi =0 (w)(n*{w) + N(w) - 2\N(w)\ 2 cos f3^j sin/3, (36) 



dp 

we find conditions of extremal values of the function Wi=\ {w,(3): 

sin(/3) = 0; 

N*(w)+N{w) Re{N(w)) (37) 

cos(P) 2|A^(w)| 2 Re(N(w)) 2 +Im(N(w)) 2 ' 

Calculations for 210 Po for the given potential have shown, that the second condition in l|37|) in a range w = 50-800 
keV is not carried out. One can explain this by that the integral J c (l, w) is smaller than the integral J c (0, w) in 10 2 -10 4 
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100 200 300 400 500 600 700 800 



Fig. 1. The Br spectra in a-decay of 210 Po: 1 is a curve, extracted from [5J by the instant accelerated model; 2 is the experimental 
data |1I4| ; 3 is the experimental data [Sj; 4 is a curve of the Br probability Wi = i in the first approximation I = 1 by our approach 
at /? = 90°; 5 is a curve of the Br probability component W^ l=1 ^ in the first approximation I = 1 at the angle 45° by our approach; 
6 is a curve calculated by us with radial integral (6) and formula dipi(r)/dr — —ipi(r)/w dV(r)/dr (it coincides with radial 
integral (7) in [5j with factor —1/w) and further normalization at the third point of data III II : 7 is a curve calculated on the 
basis of the potential (I39H - 14H 

times (that is in agreement with a condition of the convergence of the Br spectra at increasing of I). From the first 
condition in (|37|l we obtain such extremal values for the angle (3: 

/3 = 0,tt. (38) 

The Br probability with such angle values has the maximal and minimal values, and between them it varies monotonously 
for any energy of the photon emitted in the range w = 50-800 keV. One can see this also from the Table 1. 



11 Inclusion of Woods-Saxon potential into the model 

Now let's analyze, how the Br spectrum in the a-decay of the studied nucleus 210 Po is changed, if in the approach 
for calculations of the Br spectra from the potential of interaction between the a-particle and the daughter nucleus, 
pointed out in Sec. El and having the simplified barrier, to pass to a potential with the barrier, constructed on the 
basis of the account of realistic nuclear forces of interaction between the a-particlc and the daughter nucleus which is 
used in realistic nuclear models. 

With such a purpose we shall take a potential, proposed in [2S1 for description of the a-decay and synthesis of 
nuclei. Among extensive set of literature giving us different types of the a-nucleus potentials, we have given preference 
to such a paper because of here we see universal and clear approach for calculation of parameters of the potential after 
choice of needed such nucleus. In result, we suppose to obtain universal recipe for calculation of the Br spectra after 
choosing the nucleus. 

So, according to [23 (see (6)-(10)), we use such potential of interaction: 

V(r, 9,l,Q) = V c {r, 9) + V N (r, 9,Q)+Vi(r), (39) 

where 

Vc{r,6) =^(l + fj/?r 2o (0) 

V N {rAQ)= ViAZ ; Q) ro i (40) 
1 + cxp — - — 



Vi(r) 



1(1 + 1 
2mr 2 
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At current stage (with a purpose to simplify the numerical calculations of the Br spectra), for determination of the 
component Vc(r, 9) we use the formula (7) in [23 on whole region of r (without use of (8) in |25j). 
According to (14)-(20) in [2H], we calculate the parameters as follows: 

v(A, Z, Q) = -(30.275 - 0A5838Z/A 1 / 3 + 58.270/ - 0.24244Q), 
r m = 1.5268 + R, 

R = Rp(l + 3.0909/.R 2 ) +0.12434, ( , 

R p = 1.24A 1 / 3 (1 + 1.646/A- 0.1917), { ' 

t = I -QAA/(A+200), 
d = 0.49290, 



According to we see, that the parameter (3 for 210 Po is very small, that points out to a high degree of sphericity 
of this nucleus. Therefore, for calculation of the Br spectra we note the following: 

— In definition of r m (9) and R(0) we do not use (21) and (22) of t.25 (using (15) from |25|). 

— The formalism for calculation of the Br spectra, presented in Sec. l3T TTUjl . is constructed on the basis of division of 
the total wave function into its radial and angular components, i. e. in the assumption of spherical symmetry of 
the decaying nucleus. Therefore, for the nucleus 210 Po our approach for calculation of the Br spectra is applicable 
with account of realistic nuclear forces also. 

Further, we calculate the radial wave function of the decaying system for the potential (|3?)|) — (|41|) . This gives us 
the general solution for the wave function in dependence on the selected energy level for a-decay. To achieve, that 
the found solutions describe the states of the decaying system before and after the spontaneous photon emission, we 
should take into account boundary conditions in initial and final states. Here, we use such conditions: 

the initial i-state: Xi( r ~ * +oo) — > G(r) + iF(f), 

the final /-state: Xf(r = 0) = 0, { ' 

where <Pij(r) = — , F and G are Coulomb functions. 

One note, that in against of a scattering of the a-particlc on the nucleus, where as the boundary condition for 
the initial i-state a finiteness of the radial wave function <fi(r) should be used at point r = (xi( r — 0) = 0), for 
decay we choose a natural requirement, that the radial wave function tend to a spherical divergent wave in asymptotic 
region (xi{ r ~~ > 00 ) tend to a plane divergent wave). This condition gives us inevitably the divergence of the total radial 
wave function in the initial i-state at point r — (which real and imaginary parts consist from regular and singular 
solutions)! One can make sure in this by requiring a fulfillment of the continuity condition for the radial wave function 
on the whole region on its definition; or by requiring the constancy of the radial flux density, which is distinct from 
zero and is directed outside in the asymptotic region (and, therefore, it should be not zero near point r = 0, that is 
impossible to execute with null wave function at any chosen point of r). This peculiarity complicates essentially the 
calculations of the Br spectra for the a-decay in comparison with the problems of the calculation of the Br spectra 
for the scattering of charged particles on nuclei (where the essential progress has been achieved and a lot of papers 
are published). 

It turns out, that real and imaginary parts of a sub-integral function of |17| ) for the calculation of the matrix 
elements, constructed on the basis of the found solutions for the wave function for the initial and final states, tend to 
zero at point r — ► 0/ This interesting peculiarity provides the convergence of the matrix elements near point r = and, 
therefore, in the whole region of r (in the asymptotic region the convergence of the wave function is determined by the 
convergence of the Coulomb functions, considered above). Thus, we resolve the divergence problem in the calculations 
of the Br spectra in the a-decay of the nucleus 210 Po. 

The Br spectrum for the nucleus 210 Po with the a-nucleus potential i|39|) (|41|l by our approach and calculations 
is shown by the curve with number 7 in Fig. ^ From here one can see, that the new curve 7 of the Br probability is 
located very close to the curve 4 for the Br probability with the potential from Sec. with the simplified barrier (and 
also close to the curve 6 by the approaches 0E]). 



Conclusions: 

— The account of nuclear forces, essentially changed a shape of the barrier in its internal part (and essentially changed 
Br from this internal region), changes very slowly the spectrum of the total Br in the a-decay of the nucleus 210 Po 
(in comparison with the earlier obtained Br spectrum on the basis of the potential from Sec. [U| with the simplified 
barrier) . 

— This point confirms the result (obtained early on the basis of the a-nucleus potential pointed out in Sec. [5] with 
the barrier of the simplified shape) that the Br emission from the internal region till point r for barrier maximum 
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gives the very small contribution into the total Br spectrum. This conclusion coincides logically with a property 
(found on the basis of microscopic models of the nuclei with their a-decay) of leaving of the a-particle from the 
nuclear surface during the first decay stage. 

12 Conclusions and perspectives 

We present the new method of calculation of the Br spectra in the a-decay, where the angle between the directions of 
the a-particle motion (with tunneling) and the photon emission (used in the experiments [HE]) is taken into account. 
Using it, the angular spectra for the nucleus 210 Po are obtained for the first time. Now let us formulate the main 
conclusions and perspectives: 

— The method gives such a dependence of the bremsstrahlung spectrum in the a-decay of 210 Po on the angle (this 
has been obtained for the first time): 

— the first approximation at I = gives independence of the spectrum on the angle; 

— the second approximation at I = 1 gives a slow monotonous variation of the slope of the spectrum curve with 
changing the angle and without a visible change of the shape of the spectrum curve (i. e. without the appearance 
of humps and holes in the spectrum); 

— for arbitrary energy of the photon emitted in the range of w = 50-750 keV the bremsstrahlung probability is 
maximal at the angle 0° and is minimal at the angle 180°, between these angular values the bremsstrahlung 
probability varies monotonously. 

— Results for 210 Po have obtained on the basis of these approximations: 

— the bremsstrahlung process does not depend on the direction of the leaving the a-particle relatively to the shape 
of the daughter nucleus before the photon emission (this supposition has been fulfilled for 210 Po, because, in 
accordance with (see Fig. 5 on p. 33), coefficients /3° of the shape deformation for this nucleus at A = 2, 4, 6, 8 
are extremely close to zero in comparison with other nuclei with other numbers of protons and neutrons, i. e. 
210 Po is one of the most spherical nuclei); 

— the photon emission does not change the direction of the a-particle propagation (this supposition is suitable 
for the low-energy photons, and, therefore, it is applicable for analysis of all existing experimental data of the 
bremsstrahlung spectra, where one can select a region with smaller photons energies for increasing accuracy); 

— the bremsstrahlung spectra have been calculated by means of the a-nucleus potential with the simplified barrier 
pointed out in Sec. El (they coincide with the a-nucleus potential in |21| . and also with the potential with the 
external Coulomb field in 9 and in [xj) and with use of the a-nucleus potential with the barrier, pointed out 
in Sec. ^2 (see j2HJ) and constructed on the basis of realistic nuclear forces. 

— Taking into account nuclear forces in the method gives the following: 

— it does not change dependences on the angle of bremsstrahlung probability in the first and second approxima- 
tions; 

— it essentially changes the shape of the barrier in its internal region (sufficiently changes the Br from such internal 
region) and changes very little the spectrum of the total Br in the a-decay for the nucleus 210 Po at selected 
angle value. 

From here a question naturally arises: which improvement should be made in the method that results in enough 
visible changes of the Br spectrum curve, to achieve a better description of the experimental data? Note the following: 

— According to our calculations, in consideration of the possibility of the a-particle leaving at the energy of the exited 
state of the decaying system, the angle of the slope of the Br spectrum curve increases (monotonously). Apparently, 
it allows to displace the calculated Br curve (for example, by our method) essentially closer to the experimental 
data p. 

— For obtaining reliable values of the Br spectra for the a-decay one needs to achieve in the calculations the con- 
vergence of integrals for such spectra. This leads to the necessity to consider wave functions inside a large space 
region with the external boundary far enough from the nucleus. From here a new question naturally arises as to 
taking into account electrons shells of the atom with such nuclear a-decay in the calculation of the Br spectra 
in the a-decay and one can formulate the following hypothesis about the visible influence of the electrons shells 
of the atom on the total Br spectrum in the a-decay. One can note that an essential progress has been made by 
M. Amusia in the study of the Br in atomic physics |27U28| . Moreover, according to (see p. 20-21), there is 
an inevitable influence of the a-decay process at its starting time stage on the electrons shells of the atom whose 
nucleus decays. So, the a-particle during its leaving (with tunneling) deforms and polarizes these electrons shells. 
In one's turn, the changed electrons shells can correct our understanding of the real a-nucleus potential in the 
model, which should be used for the calculations of the Br spectra (one should note that these effects are still not 
studied in details). Note that these effects (partially) take place in the same space region, where it is necessary to 
use the wave functions for the calculation of the matrix elements to achieve converging values of the Br spectra. 
Therefore, it is desirable to study these effects to obtain for more accurately the total Br spectra of the a-decay. 
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— The inclusion of the a-nucleus potential from |25| in our method, deforming the decay barrier, does not essentially 
displace a point, where the a-particle starts to tunnel through the barrier. It turns out that the displacement of 
this point is much smaller in comparison with the tunneling region and even with a "mixed region", introduced in 
jSj. Therefore, after taking into account the realistic nuclear forces in the method, the interest to analyze the Br 
from these regions (with detailed study of tunneling) remains. 

— We assume that further development of the time formalism for the description of the Br in the a-decay at its 
first stage will give new abilities in the accurate description of the Br. One cannot exclude the assumption about 
appearance of "holes" in the Br spectra (see (2K])j that can allow to better describe the experimental data [2] 
3 . However, in such a case it is not clear how to connect this with the available experimental data pQ without of 
"holes" in the logical basis of our method. 

Now let's formulate the conclusions, which have a physical sense and on the basis of the calculations for 210 Po by 
our model: 

— The bremsstrahlung in the a-decay of the spherical nuclei depends on the angle extremely weakly. In taking into 
account of the nuclear forces, such dependence is not changed visibly. 

— It is not enough to take into account only one angle for the explanation of the difference between the experimental 
spectra pQ and 0EI for 210 Po (which equals to 90° and 25°, correspondingly) on the basis of our model and for 
the explanation of the difference between these experimental spectra and the calculated curves averaged by angle 
values in the approaches [312111] (that can be supposed in 0G3). 

— The small visible change of the Br spectra after taking into account of the realistic nuclear forces in our method 
confirms the result (obtained on the basis of the a-nucleus potential from Sec. with the barrier of the simplified 
shape, also see \21\ ) that the Br from the internal region till point r for barrier maximum gives the very small 
contribution into the total Br spectrum. This conclusion become natural if to take into account such a property 
(found on the basis of microscopic models of the a-decay) as the a-decay starts from leaving of the a-particle from 
the nuclear surface. 

In closing, supposing that the Br spectra in the a-decay must to change essentially at change of the angle value, 
we note, how this point can be explained, analyzing this question in theoretical and experimental aspects: 

— One can explain such angular change of the bremsstrahlung spectrum so: 

— In the a-decay of the (initially) spherical nuclei — by strong angular deformation of the decay barrier and 
continuous redistribution of the electromagnetic charge (or "nuclear polarization" likely the polarization of the 
electrons shells during tunneling of the a-particle, according to [2H] (see p. 20-21)). One can suppose, that here 
non-central forces between the a-particle and nucleons of the daughter nucleus play essential role, which exist 
in the barrier region mainly. Note, that a serious progress was achieved early in the microscopic approach in 
study of the bremsstrahlung in scattering of the nucleons and the a-particles on the light nuclei (see [30113111321 
EH1E]), m study of the bremsstrahlung in collisions between heavy ions and the nuclei (see [HE]), and in the 
non-microscopic approaches in study of bremsstrahlung induced by protons during their collisions on heavier 
nuclei (see [HE]). 

— In the a-decay of the deformed nuclei — by essential appearance of the angular anisotropy of the a-nucleus 
potential. Then one can extract an information about the shape of the nucleus from the angular bremsstrahlung 
spectra. 

— Experimental confirmation of the change the Br spectrum in the a-decay of the spherical nuclei at change of the 
angle value gives the following: 

— It will prove an existence of essential microscopic forces between the a-particle and the nucleons of the daughter 
nucleus, reinforcing the angular deformation of the barrier. 

— It will prove a visible influence of the bremsstrahlung on dynamics of the a-decay. In accordance with our 
model, it will be an experimental confirmation of the effect of variation of the barrier penetrability in result of 
the emission during tunneling of a charged particle, proposed in |20j . 

The angular analysis of the bremsstrahlung spectra gives the new additional information about the a-decay. 
Therefore, further angular experimental measurements of the bremsstrahlung spectra will have the interest. 
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